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INTRODUCTION 


'Die  Rosenthal  thermal  transient  test  for  bridgewire-type  electro  explosive 
devices  (EED)  is  used  for  nondestructive  testing  of  the  interface  between  the 
wire  and  the  explosive  (ref  1). 

The  thermal  transient  test  is  based  on  the  application  of  a  current  to 
heat  the  bridgewire,  and  then  monitoring  the  change  in  resistance  of  the  wire 
as  it  heats  to  determine  its  temperature  response.  For  this  reason,  only 
bridgewires  with  a  measurable  thermal  coefficient  of  resistance  can  be  tested. 
This  work  describes  a  simplified  testing  system  and  an  improved  method  of  data 
analysis. 


THEORETICAL  BACKGROUND 


According  to  Rosenthal's  differential  equation  for  bridgewire  temperature 
rise  as  a  lumped  thermal  system  (ref  1),  bridgewire  heat  balance  is: 

Cp  It  +  Y<t>  =  p(t)  “  l2V1+a<t>)  (D 

where 

C  =  heat  capacity,  watt-second/degree 
$  =»  temperature  increase,  °C 
t  =*  time,  second 

y  =  heat  loss  factor  (linear  thermal  conductance  term),  watt/degree 
p  =  power  input,  watts 
I  =  excitation  current,  amps 
RQ  =  initial  resistance,  ohms 

a  ■  temperature  coefficient  of  resistivity  of  bridgewire 
The  solution  for  equation  1  is: 

[1-exp  (-A] 

<D(t) - - 2 -  (2) 

2 

where  y"  =  y-I  RQa 

The  voltage  drop  v(t)  across  the  bridgewire  is: 

v(t)  =  IR(t)  (3) 

where  R(t)  =  RQ[1  +  a <}>( t )  ] . 
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By  combining  equations  2  and  3,  the  result  is: 


v(t)  =  1R  {1  +— ^  [1-exp  (-  -£-£■)]}.  (4) 

Y  P 

Ihe  bridgewire  thermal  response  from  equation  4  is  represented  in  figure  1. 

The  apparent  time  constant  is  defined  as 
C 

.  (5) 

Y 

This  time  constant  is  derived  from  the  current  wave  form,  which  is  intrinsic  to 
an  exponential  rise  in  the  bridgewire  thermal  response.  From  the  above  equa¬ 
tions  and  figure  1,  useful  thermal  parameters,  which  can  define  thermal 
behavior  of  the  bridgewire-explosive  interface,  are  derived  as  follows  (refs  2 
through  9). 

Temperature  rise  (<J>)  is  computed  from 

♦  -ifV  <» 

o 

Thermal  conductance  (y"),  modified  heat  loss  factor  due  to  feedback,  from 
al3R>2 

y'  =  — - - ,  watts/ohm  change  .  (7) 


Heat  capacity  of  the  system  (Cp)  is 

3  2 

arR 

c  =  — o~ ->  where  S  =  -r-  . 
p  S  dt 


Thermal  time  constant,  (x')  from 
t'  CD 

x'  =  -n  -/it  =  — r  =  time  to  0.63  (Vmax),  seconds 
U.oV  y 


where  t'  *=  time  to  half  voltage,  second. 
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SIMPLIFIED  TEST  SYSTEM 


A  simplified  thermal  transient  tester  was  developed  which  can  be  used  to 
control  initial  current  and  null  out  initial  resistance  of  EED.  Additional 
instrumentation  included  a  square  wave  generator,  a  digital  voltmeter,  and  two 
oscilloscopes.  One  was  used  to  monitor  the  wave  form  and  the  other  for 
recording  a  digitized  wave  form.  A  HP85  desktop  computer  was  used  to  calculate 
results  from  the  digital  recording  and  plot  data  output.  The  system  block  dia¬ 
gram  is  sL.jwn  in  figure  2. 

The  thermal  transient  test  circuit,  which  is  the  heart  of  the  test  system, 
can  be  easily  and  economically  made,  as  compared  to  commercial  instrumentation 
(for  example,  the  Model  605B  thermal  transient  test  set,  made  by  Pasedena  Sci¬ 
entific  Industries).  On  the  other  hand,  the  simplified  circuit  is  slower  and 
requires  a  skilled  operator.  Results  obtained  are  the  same  as  those  for  the 
more  complicated  automatic  system.  Either  system  is  helpful  in  understanding 
thermal  behavior  of  bridgewire-explosive  interfaces. 

The  simplified  circuit  diagram  (fig.  3)  is  basically  the  same  as  the 
standard  circuit  shown  in  figure  4.  Resistance  (R  )  and  capacitance  (C^  are 
eliminated  and  a  separate  square  wave  generator  Is  used.  By  changing  the 
grounding  location,  use  of  a  differential  amplifier  was  avoided. 

Because  it  does  not  have  R2  and  C2  in  figure  4,  the  simplified  circuit 
produces  the  various  types  of  floating  exponential  thermal  wave  forms  shown  in 
figure  5.  The  operator  must  adjust  zero  axis  and  exponential  wave  forms  to  set 
R  properly.  The  exponential  wave  form  can  be  dropped  below  the  zero  axis 
[j:ig.  5  (B)]  by  increasing  R.^  and  above  the  zero  axis  by  decreasing  R^  [fig.  5 
(A)].  When  the  values  of 1 R  and  ^  (EED  net  resistance)  are  equal,  the 
exponential  wave  form  starts  from  zero  axis  [fig.  5  ( C)  ] .  The  net  positive 
exponential  wave  form  is  proportional  to  the  value  of  AR/R0  for  the  electro 
explosive  device  and  represents  the  increasing  bridgewire  temperature.  The 
peak  value  can  be  easily  found  by  adjusting  the  wave  form  to  that  shown  in 
figure  5  (D).  The  following  operational  steps  are  used: 


1.  The  electro  explosive  device  is  mounted  in  a  shielded  test  chamber  and 
electronically  connected  to  the  shorted  test  leads. 

2.  The  square  wave  generator  is  turned  on  and  set  at  10  cps. 

3.  The  test  current  is  adjusted  by  using  the  by-pass  button  and  current 
control  knob.  It  should  be  noted  that  200  nA  for  high  sensitivity  EED  cannot 
be  used. 


4.  The  exponential  wave  form  is  displayed  on  both  oscilloscope  screens. 
Again,  it  should  be  noted  that  the  amplitude  of  the  Nicolet  digital  oscillo¬ 
scope  is  50  times  that  of  the  Tektronix  7904  oscilloscope,  since  its  input  is 
taken  from  the  output  of  the  7A13  amplifier  in  the  Tektronix  7904. 


5.  Read  hot  resistance,  equal  to  (Ri)  when  the  wave  form  corresponds  to 
figure  5  (D). 

6.  Reset  both  oscilloscopes  wave  forms  back  to  that  shown  in  figure  5 

(C). 

7.  Read  offset,  zero  time,  and  volts  needed  to  enter  data  into  computer 

8.  Read  the  voltage  and  time  when  the  cursor  on  the  oscillosope  is  placed 
on  Cj  of  the  wave-form  (fig.  6)  to  calculate  Vmax  and  stop  time,  and  enter  data 
into  computer. 

9.  Steps  7  and  8  are  used  for  thermal  time  constant,  Av  calculations,  and 
for  plotting  wave  forms.  A  typical  data  output  is  shown  in  appendix  A. 

The  original  wave  form  was  recorded  on  a  Nicolet  Explorer  III  oscilloscope 
and  stored  on  mini  floppy  disc.  Calculations  were  made  on  a  HP35  calculator. 
All  calculated  data  were  plotted  by  a  HP  7215A  plotter. 


DATA  ANALYSIS 


The  data  analysis  is  programmed  in  BASIC.  This  program  provides  calcula¬ 
tions  of  equations  6  through  9,  plotting  of  a  smoothed  thermal  exponential  wave 
form,  calculation  of  %,  and  print  out  by  paper  or  plotter.  Basically,  there 
are  seven  steps  in  the  program. 

Step  1.  Read  input  data,  which  is  needed  for  calculations  using  the  cur¬ 
sor  on  the  oscilloscope  -  tQ,  VQ,  I,  a,  R^,  RQ.  Instructions  are  included  in 
the  program. 

Step  2.  Input  data  for  plotter  -  XI,  X2,  X3,  Yl,  Y2,  Y3. 

Step  3.  Input  digital  data  which  is  recorded  on  floppy  disc  as  a  function 
of  volts  and  time  (V,t),  smooth  data,  and  plot  at  the  same  time. 

Step  4.  Calculate  shape  of  the  thermal  exponential  wave  form  to  define  S 
in  equation  8  by  moving  average  method. 

Step  5.  Calculate  Vinax  of  the  thermal  exponential  wave  form  to  define 
equations  6  and  7  by  mean  value  determination. 

Step  6.  Determine  thermal  time  constant  by  time  interpolated  for  0.63%  of 
Vmax  to  define  equation  9. 

Step  7.  Print  out  all  results.  Plot  of  %  Vmax  versus  time,  Vmax,  I,  RQ, 
a,  t,  slope,  4>,  \.  and  Cp. 

A  listing  of  the  HP85  BASIC  program  is  given  in  appendix  B.  Sample  re¬ 
sults  for  two  MK1  squib  headers  (table  A-l)  having  a  I  mil  nichrome  bridgewire, 
before  and  after  loading  with  borax  at  10,000  psi,  are  shown  in  appendix  A. 
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CONCLUSIONS 


A  simplified  thermal  transient  testing  system  has  been  developed  and 
demonstrated  for  use  in  measuring  the  thermal  coupling  between  the  bridgewire 
and  the  explosive  surrounding  it  in  EED.  It  can  be  easily  used  to  measure 
thermal  properties  and  control  the  quality  of  loaded  EED  nondestructively . 
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Figure  2.  Block  diagram  of  thermal  transient  test  instrumentation 
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TINE,  tilliMo 

TEST  #4,  PEAK  VOLTAGE  =  1.  664  millivolts 


Figure  A-4.  Test  number  4  on  MKl  squibs,  loaded  with  borax 


1©  REM  CODE  "K3MIPE" 

20  DISP  "ENTER  SAMPLE  ID" 

36  INPUT  1$ 

40  PRINT  "  " ;I$ 

50  PRINT 

60  PLOTTER  IS  705 
70  LOCATE  36,110,20,30 
8©  DIM  X$C883, Y$C80], T$C803, S$C 
1803 
90  DISP 

100  DISP  "ENTER  START  OF  X  AXIS" 
110  INPUT  Cl 
120  DISP 
130  DISP 

140  DISP  "ENTER  END  OF  X  AXIS" 
150  INPUT  C2 
160  DISP 
170  DISP 

180  DISP  "ENTER  TIC  MARKS  FOR  X 
AXIS" 

190  INPUT  C3 
200  DISP 
210  DISP 

220  DISP  "ENTER  START  OF  Y  AXIS- 
230  INPUT  C4 
240  DISP 
250  DISP 

260  DISP  “ENTER  END  OF  Y  AXIS" 
270  INPUT  C5 
280  DISF 
290  DISP 

300  DISP  "ENTER  TIC  MARKS  FOR  Y 
AXIS" 

310  INPUT  C6 
320  X1=C1 
330  X2=C2 
340  X3=C3 
358  Y 1 =C4 
360  Y2=C5 
370  Y3=C6 

380  SCALE  X1,X2,Y1,Y2 
390  T$="TEST  #3,  PEAK  VOLTAGE  = 
2.62  millivolts" 

480  S$="" 

410  X$=" TIME ,  mill isec" 

420  Y$=“PERCENT  OF  PEAK  VOLTAGE" 
430  FRAME 
440  FRAME 
450  FRAME 

460  MOVE  X1  +  CX2-XD/2  ,  Y2 
470  SETGU 

480  I PLOT  0,-70, -2 

490  SETUU 

580  LORG  6 

510  CSI2E  4,  .5 

520  LABEL  USING  “K"  ;  T$ 

538  MOVE  X1+CX2-X1 )/2, Y2 

540  SETGU 

550  I PLOT  0,2, -2 

560  SETU'J 


570  C‘  IZE  3,  .35 

580  LABEL  USING  "K"  ;  S$ 

598  AXES  X3  -'2 ,  Y3  •  X 1  .  Y 1 

688  CS { 2E  3,. 35 

610  FOR  X-X 1  TO  X2  STEP  X3 

628  MOVE  X , Y 1 

630  SE7GU 

648  IPc.OT  8,-2, -2 

650  SETUU 

668  LABLL  USING  "K"  ;  X 
670  NEXi  X 

680  MOVE  X1+CX2-X1 >/2, Y1 
690  SETGU 
700  IPLOT  0,-6, -2 
710  LORC.  4 

729  LABEL  USING  HKn  ;  X* 

73G  SETUU 
740  LORG  8 

750  FOR  Y=Y1  TO  Y2  STEP  2*Y3 
760  MOVE  XI, Y 

770  LABEL  USING  “K"  ;  Y; “  " 

780  NEXT  Y 

798  MOVE  XI,  Y1  +  <Y2-Y1  ),'2 

800  SETGU 

810  IPLOT  -6, 0,-2 

820  SETUU 

830  DEG 

840  LOIR  90 

850  LORG  6 

860  LABEL  USING  "K"  ;  Y* 

878  LOIR  O 

880  DISP  "ENTER  OFFSET  VOLTAGE  < 
millivolt)" 

890  INPUT  G 

900  DISP  “  ZERO  TIME  OFFSET  (mil 
1 isec) ” 

910  INPUT  0 

920  DISP  “  ENTER  MAX  VOLTS  IN  mi 
1 1 ivol ts" 

930  DISP  "UNCORRECTED  VALUE" 

940  INPUT  Ml 

950  DISP  "ENTER  CURRENT  (amperes 

)" 

960  INPUT  10 

970  DISP  "ENTER  COLD  RESISTANCE" 
980  INPUT  R0 

990  D I  P  "  ENTER  TEMP.  COEFF .  OF 
RESISTANCE" 

1080  INPUT  R1 
1810  Al=0*48+i 
1028  A8=A1 
1030  A2=A1+101 
1048  R3= 1 
1850  A5=4 

I960  DIM  B$C280], A$C82003 
1070  OUTPUT  702  ,:"N0" 

1888  IOBUFFER  B$ 

1098  TRANSPFR  701  TO  B*  FHS  j  CO 
UNT  230 
1100  CLEAR  7 
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1110  U0=VRL<B$E5j 83 > 

1120  H8=VRL  <  B$C  10.*  1 33  > 

1130  V1=URL(B$C14>  293 > 

1140  H1=VAL<B$C21 > 273) 

1150  OUTPUT  702  "0302“ 

1160  IOBUFFER  A$ 

1170  TRANSFER  701  TO  fi*  FHS  ;  CO 
UNT  8192 
1180  CLEAR  7 
1190  F=0 
1200  X5=0 
1210  Y5=0 
1220  N=0 

1230  FOR  I =R 1  TO  R2  STEP  2*R3 
1240  P=256*NUM<R*CI3)+NUM<R$CI+1 
3>-65536*<NUM( A$EI3>>=128) 
1250  J=I^2- . 5 
1260  U=P*V1 
1270  U=V-G^1000 
1280  U=V*1000 
1290  T=< J-H0)*H1 
1300  T=T*1000 
1310  T=T-0 
1320  X4=T 
1330  Y4=V*100/M1 
1340  X4=X4+X5 
1350  X5=X4 
1368  Y4=Y4+Y5 
1370  Y5=Y4 
1380  N=N+1 
1390  NEXT  I 
14O0  PLOT  X6,Y6 
1410  X6=X4^N 
1420  Y6=Y4/N 
1430  PLOT  X6 i Y6 
l',40  PENUP 
1450  Y6=Y6+ . 005 
1468  Y6=INT(Y6*100>/'180 
1470  OISP  Y6,X6 
1488  IF  F=1  THEN  1500 
1490  IF  Y6>63  THEN  1710 
1500  Rl=fll+R5 
1510  R2=A2+A5 
1520  Y7=Y6 
1538  X7=X6 

1540  IF  R2<R0+40*10  THEN  1200 
1550  M1=M1* .02/1000 
1560  S=f11/Tl 
1570  T2=Ml/<  I0TR0TR1  ,• 

1580  C=R1*I0'-'3*R0A2-'S 
1598  C1=R1*R0~2*I0~3/M1 
1600  PRINT  "VOLT  (max)  =  ".;M1 
1610  PRINT  “CURRENT (ai»PS>  ="  .;  18 
1620  PRINT  "R  (ohms)  =“ R0 
1638  PRINT  "COEFF  of  RES  =";P1 
1648  PRINT 

1650  PRINT  "TIME  CONST  =  ".;T1 
1666  PRINT  “SLOPECVmBx't >=";S 
1678  PRINT  "TEMP  RISE*  O'*  =".;T2 
1688  PRINT  "HEAT  LOSS  =";C1 


1690  PRINT  “HEAT  CAPACITY  ="iC 
1700  GOTO  1748 

1710  Tl=<Y6-63)/<Y6-Y7)*<X6-X7>+ 
X7 

1720  F= 1 
1730  GOTO  1200 

1740  FOR  Z=0  TO  .75  STEP  .05 

1750  L=63 

1760  PLOT  Z j L 

1770  NEXT  Z 

1780  PENUP 

1790  FOR  Z=T1- . 5  TO  T1+.5  STEP  . 
05 

1800  L=63 
1810  PLOT  Z  > L 
1820  NEXT  Z 
1838  PENUP 

1840  FOR  L=0  TO  18  STEP  .1 

1850  2=T 1 

I860  PLOT  Z,L 

1878  NEXT  L 

1880  PENUP 

1890  FOR  L=58  TO  68  STEP  . 1 

1908  Z=T 1 

1910  PLOT  Z,L 

1920  NEXT  L 

1930  PENUP 

1940  END 
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